First Observation of the Greisen-Zatsepin-Kuzmin Suppression 
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The High Resolution Fly's Eye (HiRes) experiment has observed the Greisen-Zatsepin-Kuzmin 
suppression (called the GZK cutoff) with a statistical significance of five standard deviations. HiRes' 
measurement of the flux of ultrahigh energy (UHE) cosmic rays shows a sharp suppression at an 
energy of 6 x 10 19 eV, consistent with the expected cutoff energy. We observe the "ankle" of the 
cosmic-ray energy spectrum as well, at an energy of 4 x 10 18 eV. We describe the experiment, data 
collection, analysis, and estimate the systematic uncertainties. The results are presented and the 
calculation of the statistical significance of our observation is described. 

PACS numbers: 98.70.Sa, 95.85.Ry, 96.50.sb, 96.50.sd 



In 1966, Greisen [![, and Zatsepin and Kuzmin 0, 
proposed an upper limit to the cosmic-ray energy spec- 
trum. Their predictions were based on the assumption of 
a proton dominated extra-galactic cosmic-ray flux which 
would interact with the photons in the cosmic microwave 
background (CMB) via photo-pion production. From the 
temperature of the CMB and the mass and width of the 
A + resonance, a "GZK" threshold of ~ 6 x 10 19 eV was 
calculated, and a suppression in the cosmic-ray flux be- 
yond this energy (commonly called the GZK cutoff) was 
predicted. This is a strong energy-loss mechanism that 
limits the range of cosmic protons above this threshold 
to less than ~ 50 Mpc. 

Several earlier experiments [3, EL IS HI have reported 
the detection of one event each above 10 20 eV. A con- 
tinuing, unbroken energy spectrum beyond the predicted 
GZK threshold was later reported by a larger experiment, 
the Akeno Giant Air Shower Array (AGASA) 0,0. 

The High Resolution Fly's Eye (HiRes) experiment was 
operated on clear, moonless nights over a period of nine 
years (1997-2006). During that time, HiRes collected a 
cumulative exposure more than twice that collected by 
AGASA above the GZK threshold. The HiRes experi- 
ment observes cosmic rays by imaging the extensive air 
shower (EAS) generated by a primary cosmic ray. Ul- 
traviolet fluorescence (UV) light is emitted by nitrogen 
molecules in the wake of the EAS and collected by our 



detector. 

Forty years after its initial prediction, the GZK cutoff 
has been observed for the first time by the HiRes ex- 
periment. In this article we describe our measurement 
of the flux of cosmic rays, the resulting cosmic-ray en- 
ergy spectrum, our analysis of this spectrum to infer the 
existence of the cutoff, and our estimate of systematic 
uncertainties. 

The HiRes project has been described previously 0, 
Hot ]. The experiment consists of two detector stations 
(HiRes-I and HiResTI) located on the U.S. Army Dug- 
way Proving Ground in Utah, 12.6 km apart. Each sta- 
tion is assembled from telescope modules (22 at HiRes-I 
and 42 at HiRes-II) pointing at different parts of the sky, 
covering nearly 360° in azimuth, and 3°-17° (HiRes-I), 
and 3°-31° (Hires-II) in elevation. Each telescope mod- 
ule collects and focuses UV light from air showers using 
a spherical mirror of 3.7 m 2 effective area. A cluster of 
256 photomultiplicr tubes (PMTs) is placed at the focal 
plane of each mirror and serves as the camera for each 
telescope. The field of view of each PMT subtends a one 
degree diameter cone on the sky. 

HiRes data analysis is carried out in two ways. In 
monocular mode, events from each detector site are se- 
lected and reconstructed independently. The combined 
monocular dataset has the best statistical power and cov- 
ers the widest energy range. The dataset consisting of 
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events seen by both detectors (stereo mode data) has the 
best energy resolution, but it covers a narrower energy 
range and has less statistics [ll|. This article presents the 
monocular energy spectra from our two detectors. 

The photometric calibration of the HiRes telescopes 
has been described previously [12|. It is based on a 
portable, high-stability (~ 0.5%) Xenon flash lamp car- 
ried to each mirror on a monthly basis. Relative nightly 
calibrations were performed using Yag laser light brought 
to each cluster of PMT's through optical fibers. In ad- 
dition, the overall optical calibration of the HiRes detec- 
tors is validated by reconstructing scattered light from 
a pulsed laser fired into the atmosphere from locations 
that surround, and are within ~ 3.5 km, of the two de- 
tector sites. We achieve ~ 10% RMS accuracy in our 
photometric scale. 

We monitor the UV transmission properties of the at- 
mosphere to make a correction for the attenuation of flu- 
orescence light. Steerable lasers fire patterns of shots 
that cover the aperture of our fluorescence detectors, and 
the detectors measure the intensity of the scattered light. 
The most important parameter we measure is the verti- 
cal aerosol optical depth (VAOD). The mean value of the 
VAOD is 0.04 with an RMS variation of 0.02. An event 
at 25 km from a HiRes detector has an average aerosol 
correction of ~15% upward in energy. Because 2.5 years 
of early HiRes-I data were collected before the lasers were 
deployed, the spectra presented here are calculated using 
a constant-atmosphere assumption, using the measured 
average value for the VAOD. We have tested this assump- 
tion by calculating the energy spectrum from our later 
data, using the actual hourly measurements. Comparing 
the resulting spectra from the two methods, we obtain 
flux values that agree to within a few percent 

Another important parameter in our analysis is the flu- 
orescence yield (FY): the number of photons generated 
per ionizing particle per unit path length. FY measure- 
ments have been made by several groups 3, HI, 3, 13] • 
For the energy spectrum determination used in this pa- 
per, we have used the spectral shape of Bunne r 11411 and 
the integral yield reported by Kakimoto et al. [15j |. Our 
systematic studies have shown that this set of assump- 
tions produces absolute fluorescence flux values that are 
equal, within ~ 6%, of those obtained using a fit to all 
the results cited fl8| |. 

The details of HiRes event selection have been de- 
scribed previously An additional cut on the 
distance to showers has been applied in the HiRes-II data 
collected after those shown in [20j. This cut is applied 
to make the aperture (defined as the product of collec- 
tion area and solid angle) calculation more robust. The 
event reconstruction procedure begins with the determi- 
nation of the shower axis. A plane containing the axis of 
the shower and the detector, the shower-detector plane 
(SDP) is determined from the pointing direction of trig- 
gered PMTs. For the HiRes-II monocular dataset, the 



PMT times are then used to find the distance to the 
shower and the angle, ip, of the shower within the SDP. 
This timing fit measures ip to an accuracy of ~ 5° RMS. 

The number of shower particles as a function of atmo- 
spheric depth is then determined. This calculation uses 
the FY and corrects for atmospheric attenuation. We fit 
this shower profile to the Gaisser-Hillas function 2l| , af- 
ter having subtracted scattered Cerenkov light produced 
by the air shower particles. This profile fit yields both 
the energy of the shower and the depth at the shower 
maximum, A max . A typical HiRes profile is displayed in 
[l2T |. The energy resolution of the HiRes-II detector is 
about 12% at high energies. 

The HiRes-I detector, with its limited elevation cover- 
age, does not typically observe enough of the shower for 
a reliable timing fit. For this reason the HiRes-I monoc- 
ular reconstruction combines the timing and profile fits 
in a profile-constrained fit (PCF). The PCF reconstructs 
ip with an accuracy of ~ 7° RMS. The PCF has been 
validated by comparing the PCF energies to those found 
using stereo geometries in that subset of the data ob- 
served by both detectors as shown in Figure [TJ The en- 
ergy resolution of the HiRes-I detector is about 17% at 
high energies. 
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FIG. 1: HiRes-I energies calculated with the event geom- 
etry reconstructed in monocular mode using the profile- 
constrained fit vs. the energy reconstructed in stereo mode. 

Finally, a correction is made for the energy carried by 
shower components which do not deposit their energy 
in the atmosphere. This correction includes primarily 
the energy of neutrinos and muons that strike the earth. 
The correction is calculated using shower simulations in 
CORSIKA 22| with hadronic interaction simulated by 
QGSJet [23|]. The correction is ~ 10%. Simulations using 
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Sibyll 24] find a correction within 2% [l3j of that found 
via QGSJet. 

The measurement of the cosmic-ray flux requires a re- 
liable determination of the detector aperture. The aper- 
ture of the HiRes detectors has been calculated using a 
full Monte Carlo (MC) simulation. The MC includes sim- 
ulation of shower development (using CORSIKA), fluo- 
rescence and Cerenkov light production, transmission of 
light through the atmosphere to the detector, collection 
of light by the mirrors, and the response of the PMTs, 
electronics and trigger systems. Simulated events are 
recorded in the same format as real data and processed in 
an identical fashion. To minimize biases from resolution 
effects, MC event sets are generated using the published 
measurements of the energy spectrum [25( and composi- 
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To ensure the reliability of the aperture calculation, 
the MC simulation is validated by comparing key distri- 
butions from the analysis of MC events to those from the 
actual data. Several of these comparisons were shown in 
reference [2f|. Two comparisons are especially notewor- 
thy. The data-MC comparison of the distances to showers 
shows that the simulation accurately models the cover- 
age of the detector. The comparison of event brightness 
shows that the simulations of the optical characteristics 
of the detector, and of the trigger and atmospheric con- 
ditions, accurately reproduce the data collection environ- 
ment. The excellent agreement between the observed and 
simulated distributions shown in these cases is typical 
of MC-data comparisons of other kinematic and physical 
quantities, and this agreement demonstrates that we have 
a reliable MC simulation program and aperture calcula- 
tion. Figure^] shows the result of the aperture calculation 
for both HiRes-I and HiRes-II in monocular mode. 
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FIG. 2: The apertures (defined as the product of collection 
area and solid angle) of the HiRes-I and HiRes-II detectors 
operating in monocular mode. 

Figure [3] shows the monocular energy spectra from the 
two HiRes detectors [30| . The data included in the figure 
were collected by HiRes-I from May, 1997 to June, 2005, 
and by HiRes-II from December, 1999 to August, 2004. 
Figure [3] shows the flux multiplied by E 3 , which does not 



change the statistical interpretation of the results but 
highlights features more clearly. Two prominent features 
seen in the figure are a softening of the spectrum at the 
expected energy of the GZK threshold of 10 198 eV, and 
the dip at 10 186 eV, commonly known as the "ankle". 
Theoretical fits to the spectrum [3l| show that the ankle 
is likely caused by e + e _ pair production in the same in- 
teractions between CMB photons and cosmic-ray protons 
where pion production produces the GZK cutoff. The ob- 
servation of both features is consistent with the published 
HiRes results of a predominantly light composition above 
10 18 eV El. 
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FIG. 3: The cosmic-ray energy spectrum measured by the 
HiRes detectors operating in monocular mode. The spectrum 
of the HiRes-I and HiRes-II detectors are shown. The highest 
two energy bins for each detector are empty, with the 68% 
confidence level bounds shown. The spectrum of the AGASA 
experiment is also shown 0, H|. 

At lower energies, the cosmic-ray spectrum is well fit by 
a piece-wise power law model. A similar fit also gives an 
excellent representation of the spectrum in Figure [3] The 
three straight line segments shown represent the result 
of a fit of the measured flux to a triple-power law. The 
fit contains six free parameters: one normalization, the 
energies of two floating break points, and three power 
law indices. 

We performed a binned maximum likelihood fit [32| to 
the data from the two detectors. The fits include two 
empty bins for each monocular dataset. We found the 
two breaks at \ogE (E in eV) of 19.75 ±0.04 and 18.65 ± 
0.05, corresponding to the GZK cutoff and the ankle, 
respectively. When the datasets were made statistically 
independent by removing events seen by both detectors 
from the HiRes-I dataset, we obtained a x 2 of 35.1 in this 
fit for 35 degrees of freedom (DOF). In contrast, a fit to a 
model with only one break point, while able to locate the 
ankle (at the same energy), yielded a x 2 /DOF=63.0/37 

m. 

A measure of the significance of the break in the spec- 
tral index at 10 198 eV can be made by comparing the 
actual number of events observed above the break to the 
expected number for an unbroken spectrum. For the lat- 
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ter, we assume the power law of the middle segment to 
continue beyond the threshold. From the independant 
HiRes exposures (with events seen by both detectors re- 
moved from HiRes-I), we expect 43.2 events above 10 19 8 
eV from the extrapolation, whereas 13 events were actu- 
ally found in the data. The Poisson probability for the 
observed deficit is 7 x 10~ 8 , which corresponds to 5.3 
standard deviations. We conclude that we have observed 
the GZK cutoff with a 5 standard deviation significance. 

One question that remains is whether the sources of ex- 
tragalactic UHE cosmic rays have properties that could 
change the GZK energy. A study by V. Berezinsky et 
al. [3lJ found that the density of sources in the local 
area should change the power law of the energy spectrum 
above the GZK cutoff, but not the GZK energy itself. 
The average power law of the sources could change the 
GZK energy somewhat, but the method suggested 



by V. Berezinsky and S. Grigorieva [34j provides a test of 
whether a break is the GZK cutoff independent of power 
law over a wide range. -E1/2 refers to the energy at which 
the integral energy spectrum falls to half of what would 
be expected in the absence of the GZK cutoff. To calcu- 
late £1/2 we used the HiRes monocular energy spectra 
and the integral of the power law spectrum used above to 
estimate the number of expected events above the break. 
We find E1/2 — 10 19 - 73±0 07 . Berezinsky and Grigorieva 
predict a robust theoretical value for Ei/ 2 of 10 19 - 76 eV 
for a wide range of spectral slopes 34J . These two values 
are clearly in excellent agreement, supporting our inter- 
pretation of the break as the GZK cutoff. 

We measure the index of the power law to be 3.25±0.01 
below the ankle, 2.8f ± 0.03 between the ankle and the 
GZK cutoff, and 5.1 ± 0.7 above the GZK cutoff. 

For the monocular analyses, the main contributions to 
the systematic uncertainty in the energy scale and flux 
measurements are: PMT calibration (10%), fluorescence 
yield (6%), missing energy correction (5%), aerosol com- 
ponent of the atmospheric attenuation correction (5%), 
and mean energy loss rate estimate (the flux of fluores- 
cence photons is proportional to the mean dE/dx of the 
particles in the shower [35j) (10%). Since these uncer- 
tainties arise from very different sources, we add them 
in quadrature, giving a total energy scale uncertainty of 
17%, and a systematic uncertainty in the flux of 30%. 

In summary, we have measured the flux of ultrahigh 
energy cosmic rays with the fluorescence technique, in 
the energy range 10 17 2 to above 10 20 5 eV. We observe 
two breaks in the energy spectrum consistent with the 
GZK cutoff and the ankle. The statistical significance of 
the break identified with the GZK cutoff is 5 standard 
deviations. We measure the energy of the GZK cutoff to 
be (5.6 ± 0.5 ± 0.9) x 10 19 eV, where the first uncertainty 
is statistical and the second is systematic. 
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